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GLOBAL DYNAMICS FOR STEEP SIGMOIDAL NONLINEARITIES IN 

TWO DIMENSIONS 

TOMAS GEDEON, SHAUN MARKER, HIROSHI KOKUBU, KONSTANTIN MISCHAIKOW, 

AND HIROE OKA 


1. Introduction 


This paper discusses a novel approach to obtaining mathematically rigorous results on 
the global dynamics of ordinary differential equations. The motivation is twofold. The 
first arises from applications, in particular the study of regulatory networks. Because of 
their centrality in the study of systems biology anything beyond the most cursory com- 
ments are beyond the scope of this introduction. Instead we refer the interested reader to 
Albert et al.l 2013]. The second arises from ongoing work of the authors to develop a math¬ 
ematical framework, which we call a Database for Dynamics, that provides a computation¬ 
ally efficient and mathematically rigorous analysis of global dynamics of multiparameter 
nonlinear systems. Rea sonable success h a s been obta i ned i n the context of nonlinear sys ¬ 
tems generated by maps Arai et ^ 2009l | , Bush et al. (2012 |. Bush and Mischaikow 2014l | , 


Bush et al.l [2015l |. However, extending these methods to ordinary differential equations is 


proving to be technically challenging since the rigorous evaluation of a map needs to be 


replac e d by the rigo r ous ev aluation of solutions to a differential equation Mischaikow et al 
2 OI 5 I ]. Mivaii et al.l 2015j |. We return to this topic in the conclusion. 


A regulatory network can be represented as an annotated directed graph. The vertices 
represent a regulatory object, e.g. a protein, and an edge from node m to node n indicates 
that m directly regulates n. The annotation on this edge indicates whether m activates (up 
regulates) or represses (down regulates) n. It is natural to think of modeling the dynamics 
of this system via an ordinary differential equation. In the context of gene regulatory 
networks, the associated proteins have a natural decay rate, and thus given a network with 
N genes one is led to a system ordinary differential equation of the form 

(1) Xn = -'lnXn +fn{x)., n=l,...,iV 

where 7 ^ > 0 and fn is independent of Xm if there is no edge from node m to node n. In 
the context of applications at this level of generality there is little that one can hope to say. 
However, it is typically assumed that the interactions have switch like behavior and thus 
the nonlinearities used to model interactions between individual nodes are often assumed 
to have a sigmoidal shape. A common construction is based on Hill functions 

( 2 ) 


e^ + : 


or 
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where the former and latter expressions are used to model activation and repression, re¬ 
spectively. However, it is important to keep in mind that ther e are other models that are 
probably more representative of the underlying biochemistry Bintu et al. 2005 Even 
for N of moderate size the analysis of ([T]) with arbitrary Hill functions is intractable. A 
standard simplification is to let A: —)• oo in which case one obtains nonlinearities that take 
the form of piecewise constant functions. We refer to this as a switching system (for a 
precise definition see Section [3]) and denote it by 

(3) i; = —Ex-|-A(x), x E 

We view ([3]) as a computational model, the purpose of which is to give us insight into 
the behavior of the biologically motivated model ([l]) in which the nonlinearities fn are 
Lipschitz continuous though a particular analytic form is not known. In particular, we are 
not concern ed with identifying solutions to ([3]) . 

We recall Albert et al.l 20131 ] that there is a long tradition of associating state transi¬ 
tion diagrams, which take the form of directed graphs, to regulatory networks. The paths 
through the state transition diagram are then used to represent the dynamics of the net¬ 
work. As is described in Section [6] we use ([3]) to define a particular choice of state transition 
diagram. Because our focus is on dynamics we find it convenient to represent this state 
transition diagram as a combinatorial multivalued map V ^ V. In this notation, V 
denotes the set of vertices in the state transition diagram and there exists a directed edge 
tt — >• X in the state transition diagram if and only if x € iF{u). 

The number of elements in V can grow rapidly as a function of the size of the regulatory 
network. In particular, in the approach we take here, if 0{n) denotes the number of out 
edges at node n, then the size of V is of order 

N 


n(0(n) + l). 


n=l 

Cataloguing all the paths in a graph of this size is not practical. However, there are efficient 
(both in time and memory) graph algorithms that allow one to identify essential dynamical 
structures: the recurrent dynamics, i.e. the nontrivial strongly connected components of 
iF] and the gradient-like dy namics, i.e . the r e achability, defined by paths in iF, between the 
recurrent components (see Arai et ^ 20091 ] . Bush et al. 2012| ] and references therein and 
Berenguier et al. |2013| ] for an application of these techniques in the context of state tran¬ 
sition diagrams). We encode this information in the form of a Morse graph, MG(J^). This 
is the minimal directed acyclic graph such that each nontrivial strongly connected compo¬ 
nent is represented by a distinct node and the edges indicate the reachability information 
inherited from iF between the nodes. 

The switching system ([3]) has an N dimensional pha se space, but an i y-|-3E dimensional 


parameter space {E denotes the number of edges). In I Cummins et al.l 2015bl ] we describe 


a set of algorithms that allows us to decompose parameter space into semi-algebraic sets, 
such that on each set the state transition diagram E is constant and hence the global 
dynamics as described by the Morse graph M G(J^) is valid for eac h parameter value in the 
set. The algorithms have been implemented Harker et al. 20151 ] and thus for moderate 
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sized N we have the capability of describing via the Morse graphs the global dynamics 
associated with the state transition diagrams for all parameter values. 

As is emphasized above we believe that from the biological perspective © provides a 
more realistic model than m- Therefore, to justify the biological relevance of the combi¬ 
natorial computations described above that are based on ([3]), it is important to be able to 
demonstrate that MG(T') provides correct and meaningful information about the dynamics 
of ([H). Since a Morse graph is a directed acyclic graph it generates a poset. Our goal is to 
translate the poset structure associated with MG(T') into information about the structure 
of invariant sets for a flow generated by a smooth differeritial eq uation. For this we make 
use of Morse decompositions as defined by Conley Conlev 1978| | . Recall that given a con¬ 
tinuous flow (/?: M X X — )• Ai defined on a compact metric space X a Morse decomposition 
of X consists of a finite collection of mutually disjoint compact invariant sets called Morse 
sets indexed by a partially ordered set (P, <) with the property that if 

xEX \ IjM(p), 

pGP 

where M(p) denotes the Morse set indexed by p, then 

a{x,ip) C M{p) and uj{x,(p) C M{q) 


where q < p. 

The primary goal of this paper is Theorem 18.21 which roughly states that given a reg¬ 
ulatory network we can use ([3]) to construct a state transition diagram for which we can 
efficiently compute a Morse graph MG(T') from which we can determine a Morse decompo¬ 
sition for the dynamics defined by a smooth system of the form ([T]). The following outline 
indicates the tools and constructions that are used to obtain a proof of Theorem 18.21 

We begin in Section [2] with a brief description of Conley theory. Section 12.11 contains 
elementary ideas from lattice and poset theory and a statement of Birkhoff’s theorem 
that relates finite distributive lattices and finite posets. Sections 12.21 and 12.31 presents 
the necessary definitions of Conley theory in the settings of combinatorial and continuous 
dynamics, respectively. Section 12.41 contains Theorem 12.91 which provides the theoretical 
framework by which we translate the information from the combinatorial dynamics to the 
continuous dynamics. 

As indicated above Section [3] presents the definition of a switching system ([3]) and associ¬ 
ated notation and definitions. The material in this section is restricted to two-dimensional 
systems, but can straightforwardly be extended to systems of arbitrary finite dimension. 
In Section 0] (5-constrained continuous switching systems 

(4) x = —Tx + f^^\x) 

are defined. The definition begins with a given switching system, which determines F and 
a positive number 5, which indicates the width of a collar around the lines of discontinuity 
of A. The function is obtained by replacing A by a continuous function on this collar. 
Again, this construction is done in the setting of M^, but can be extended to M”. 
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Theorem 18.21 states that the MG(J-') derived from the switching system ([3]) determines 
a Morse decomposition for any associated 5-constrained continuous switching system (jl]) 
for any 0 < (5 < 5*, where 5* is explicitly determined by T and A (I32p . The strategy of the 
proof is as follows. We use the results of Section [2] and in particular Birkhoff’s theorem 
to pass from the poset structure induced by MCPT”) to a lattice of attractors for the state 
transition graph We use Theorem 12.91 to guide the construction of a lattice of forward 
invariant sets for the 5-constrained continuous switching system. Birkhoff’s theorem is 
then once more employed to identify the poset structure of the Morse decomposition for 
(5-constrained continuous switching system. 

The state transition graph T associated with a switching system ([3]) is defined in Sec- 
tion El This constructi on is presented in the setting of but has been extended to M” 
' 2ni5bl |. 


Cummins et al 


The major technical work of this paper is to construct a lattice of trapping regions 
for the (5-constrained continuous switching system (jl]) that is isomorphic to the lattice 
of attractors for T. This is done in two steps, both of which are restricted to M^. The 
first is to construct elementary regions in that we call tiles and chips, which are used 
to construct the trapping regions. Tiles are related to the definition of and defined 
in Section U) Chips are defined in Section [5] along with a proof that if the constraints 
on (5 imposed by (l3^ are satished, then the flow of (H]) is transverse along the edges of 
interest of tiles and chips. The second step involves the construction of the desired trapping 
regions. This is done in Section [ 7 ] in slightly more generality than needed as we construct a 
trapping regions associated with any forward invariant set of T". There are two important 
remarks that need to be made about the content of Section [71 First, it is a tedious case by 
case local analysis of the constructing and verihcation of the trapping regions and hence 
too cumbersome to generalize to higher dimensions. Second, we do not know of counter 
examples to the construction for higher dimensional systems, thus we believe that with an 
alternate proof it might be possible to generalize Theorem 18.21 to n-dimensional systems. 

The formal statement and proof of Theorem 18.21 along with related results are presented 
in Section [8l 


2. Conley Theory 


Our proof that efficient graph theoretic computations can lead to rigorous mathematical 


results for smooth switching s ystems is ba s ed on new develop ments 


presented in a series of papers iKalies et al 


i ed on new aevelop mt 
(2nn,l I2ni4 l2ni5al lb|. 


We 


in Conley theory as 
review the essential 


ideas of these results along the lines in which we employ these ideas in this paper. We 
begin with a brief review of posets and lattices, describe Conley theory first in the context 
of directed graphs (combinatorial dynamics), and then in the setting of continuous flows on 
compact metric spaces. Finally, we state a theorem using this language that provides the 
framework by which the state transition graph leads to mathematically rigorous statements 
about the global dynamics for a smooth switching system. 
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2.1. Posets and Lattices. We assume the read e r is familiar with the concepts of pa rtially 


ordered sets (posets) and lattices (see Roman 2008( |. Davev and Priestley 2002l |h but 


review some fundamental concepts as a means of establishing notation. 

Given a poset (P, <), D C P is a down set of P if p G P, g G D, and p < q implies that 
p £ D. The set of down sets of P is denoted by 0(P) and is a lattice under the operations 
of union and intersection. In fact, 0 defines a contravariant functor from the category of 
posets to the category of lattices. 

Let (L, V, A,0,1) be a bounded distributive lattice where 0 and 1 denote the minimal 
and maximal elements, i.e. 


0 A U = 0 and 1 A U = U 


for all 17 G L. The lattice algebra induces a partial order on its elements as follows. Given 
U,V gL 

(5) a UAV = U then U <V. 

Recall that a nonzero element U G L is join irreducible if 7/ = AW B implies that U = A 
or U = B. We denote the set of join irreducible elements of L by J^(L). Observe that in a 
finite lattice if 7/ is a join irreducible element, then it has a unique predecessor with regard 
to the partial order ([5]). We denote this unique predecessor by 

(6) U . 

Using the partial order ([5|) (J^(L),<) is a poset and, more generally, defines a con¬ 
travariant functor from the category of bounded distributive lattices to the category of 
finite posets. 

Theorem 2.1 (Birkhoff’s theorem). Let (L,V,A) he a finite distributive lattice and let 
(P;<) be a finite poset. Then, 

(i) 0(J^(L)) is lattice isomorphic to L. 

(ii) J^(0(P)) is poset isomorphic to P. 

As is made clear below our interest in Birkhoff’s theorem lies in the fact that it guarantees 
that one can represent the same information either in a poset or a lattice and that there is 
a well defined transformation between the two. 


2.2. Combinatorial Conley Theory. As indicated in the Introduction, to emphasize 
the fact that we are interested in dynamical structures we represent a directed graph as a 
combinatorial multivalued map J-: V ^ V, where V is the finite set of vertices and there 
is a directed edge z/ —)• z/' if and only if v' G B{v). Note that we allow self edges in our 
directed graph thus it is possible that v £ We use the notation z/ ^ z/' to indicate the 

existence of a path from u to u'. Using the multivalued map notation, z^ z^' is equivalent 
to the statement that there exists re G N such that v' £ J-'”(z/). Backward paths in the 
graph can be associated with reversal of time. With this in mind define : V —)• V by 

v' £ if and only if v £ F{v'). 
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Elements E V belong to the same strongly connected path component oi J- if ly 
v' and v' ^ v. Since we allow self edges it is possible that a strongly connected path 
component consists of a single vertex with a self edge. We refer to a strongly connected 
path component of as a Morse set of J- and denote it by C V. The collection of all 
strongly connected path components of T is denoted by 

MD(A) := {M{p) C V I p€ P} 

and forms a Morse decomposition of J-. We impose a partial order on the indexing set P 
of MD(A') by defining 

q < p if there exists a path in A from an element of Ai{p) to an element of A4(q}. 

Definition 2.2. The Morse graph of A, MG(A'), is the Basse diagram of the poset (P, <). 
We refer to the elements of P as the Morse nodes of the graph. 


As discussed in the introduction given A: V ^ V identification of Morse sets and the 
Morse graph is computationally feasible. Thus, these are the objects that we extract from 
the state transition diagram. However, we know of no direct means of transferring knowl¬ 
edge of the Morse graph to an associated continuous system of the form ([3|). As indicated 
above Birkhoff’s theorem guarantees that we will not lose information by considering the 
lattice of down sets 0(P), where (P,<) is the down set that defines the Morse graph. To 
identify this lattice in the directed graph A: V ^ V recall the following concept. 


Definition 2.3. A set AA C V is forward invariant under A if A(AA) C Af. 

The collection of all forward invariant sets of A is denoted by Invset^(A) and as discussed 


m 


Kalies et al.l . l201^ . Section 2] is a bounded distributive lattice where 0 := 0 and 1 = V. 


Given Af E lnvset'*'(A) define 

Af° :={u£Af\ A-^u) n AA / 0} . 

Of central interest is the following special type of forward invariant set. 


Definition 2.4. A set A C V is an attractor for A if A(A) = A. 

Observe that given an attractor A, = A. 

The collection of all attractors in V under A is denoted by Att(A) and as discussed 
Kalies et ah . l2015al . Section 2] is a bounded distributive lattice where 0 := 0 and 


m 


1 = max {A I A E Att(A)}. Furthermore, given Aq, Ai E Att(A) the lattice operations are 
defined by 


( 7 ) 

and 


Aq V Ai :— Ao U Ai 


(8) Ao A Ai := max {A E Att(A) | A C Aq O Ai} . 

Given a Morse set A4(p) let 

4,(Af(p)) = {u E V I 3A E M.{p) such that A ^ u} . 
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Note that KM-iv)) £ Att(J^) and, in fact, {|(A^(p)) | p E P} generates Att(J^). Thus 
Att(T') = 0(P) fsee lKalies et alJ 2015b | for details). 

We make use of the following result, which follows via a finite induction argument from 
the definition of a forward invariant set. 


Proposition 2.5. Assume M E lnvsef*’(T'). If u ^ M and u 




then v' E M . 


2.3. Continuous Conley Theory. We assurne the reader is familiar with basic concepts 
from the theory of dynamical systems Robinsoi] (l999| | . Let 99 : RxX^Xbea continuous 
flow defined on a compact metric space. Let lnvset(X, ip) denote the collection of invariant 
sets in X under (/?. 

Definition 2.6. A compact set C A is an attracting neighborhood if uj{N,{p) C int(A) 
where uj{N, (p) denotes the w-limit set of N under cp. 

The se t of all attracting n eighborhoods in X under tp is denoted by ANbhd(A, ip) and as 
shown in Kalies et al. (2014 ] is a bounded distributive lattice (but in general is not finite). 
In this paper we make use of a special class of attracting neighborhoods. 

Definition 2.7. Assume the flow ip is gene rated by a di fferential equation x = f{x), 
X E Let A C be a regular closed set Walker 1974l | whose boundary is made up 
of a finite number of straight line segments. For each closed boundary edge e of A let rig 
denote the outward normal vector. We say that A is an attracting block if f{x) ■ Uf, < 0 
for every x E e. 


Recall that an invariant set A E lnvset(A, (/?) is an attractor for ip if there exists an 
attracting neighborhood such that A = uj{N,ip). The dual repeller of an attractor A is 
defined to be 

A* := {x E A I a;(x, 99) n A = 0} . 

In what follows we assume that A is a finite bounded sublattice of ANbhd(A, (^9). Recall 
that J^(A) forms a poset. Let Q —)■ J^(A) be a poset isomorphism, i.e. we are using the 
poset (Q, <) to index the elements of J^(A). Define /r: Q —)• lnvset(A, (^ 9 ) by 


p{q) := lnv{L{q),ip) n Inv L{q),ip 


As is shown in Kalies et al. (2015bj | p{Q) is a Morse decomposition of A under ip, i.e. it is 
a collection of mutually disjoint compact invariant sets with the following property: if 


geQ 

then there exists q,q' £ Q such that 

a{x,ip) C /u(g) and uj{x,ip) C p{q') 

and furthermore q' < q under the partial order on Q. The Hasse diagram of Q is the Morse 
graph associated to the Morse decomposition. 
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2.4. The Translational Theorem. 


Definition 2.8. Given a poset (P, <p) define (P, <p) to be the poset where P = P U {p} 
and <p restricted to P equals <p along with the additional relations p <p p for all p E P. 


The primary result of this paper is a corollary of the following theorem iKalies et al 


2015bt]. 


Theorem 2.9. Letip: M x X —X be a flow on a compact metric space. Let (N,A,V,0,1) 
be a finite distributive lattice that satisfies the following properties: 

N C lnvsef''(X, yj). 

( 2 ) 0 = 0 , 

(3) N ^N' d N ^^N' andNy N' = Nd N'. 

Let the poset (P,<) = (J^(N),<) be an indexing set for J^(N). For each p E P define 


i(Inv(l.v)r 


Inv ( N{p),ip 


i/p E P 
ifp = p. 


Then the collection of invariant sets M{p), p G P defines a Morse decomposition of p and 


<p is an admissible order. 


3. Two-dimensional Switching Systems 

In this section we provide a formal definition of a general two-dimensional switching sys¬ 
tem and provide elementary results about the associated dynamics. We begin with two sets 
of non-negative real numbers H := {.^j | i = 0,..., / 1} and H := {rjj | j = 0,..., J -|- 1} 

that we refer to as threshold values, with the property that 

0 = ^0 < < • • • < // < //+! = oo 

0 = rjo < m < ■ ■ ■ < hJ < hJ+i = oo. 

Let 

(9) n := {(Ci, Pj) I i = 0,..., j = 0,... , J} C [0, oo)^. 

In addition we assume that we are given partitions of the sets of threshold values 

and H = H^UH2 . 

Of primary importance is the following collection of open rectangles, called cells, defined 
in terms of the thresholds as follows: 

K, := {k(z, j) := {ii,ii+i) X C (0,oo)^ | z = 0,...j = 0,..., j} . 

Definition 3.1. The switching system S = S(r, A, H^, H^) is defined to be the 

system of differential equations 

(10) X = —Tx -b A(x), X E K C (0, oo)^ 

k£K 
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where 


and 


r = 


7i 0 
0 72 


7i > 0 


A(x) = 


Ai(x) 

A2(x) 


is constant on the cells /C. Fnrthermore, A satisfies the following constraints for all i G 

{0,..., I}, j G {0,..., J}, 


( 11 ) 


For the sake of simplicity we assume that if k* ^ Kj, then 

A(Kj) ^ A(kj). 

Fnrthermore, setting 

(12) $(ac) := F-^A(k) 

we assnme that for all ac G /C 




Ai(K(i, j)) 




Vj G Hi 



rjj G H2 


Ai(K(i, j)) 


(13) <1>i(k) 0H and <h 2 (K) 0 H. 

A consequence of assumption (jl3|) is that for a fixed switching system S 

(14) ^ = ^(E);= min {|$i(k) - |$ 2 (n) - r?|} > 0. 

kS/C, t;GH 

Two other positive constants that are used later in the paper are 

(15) p = p{T.):= inax {|4 >i(k) - ^i|, |4 >i(k) - |4>2(k) - |d'2(«:) - , 

which provides a measurement of the maximal displacement of the attracting fixed point 
^{k) from the cell k, and 

(16) 7 = 7 (S) := min | , —1 < 1. 

172 71 J 

Given a particular cell k = K{i,j) G JC there is an associated affine vector field, which 
we call the K-equation, given by 

(17) x = —rx + A(K). 

We denote the flow generated by the K-eqnation by Observe that 4 >(k), as defined in 
m, is an attracting fixed point for the K-equation. 
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K2 

'k 

Kl 

'k 

k: k: k: 
K 4 


Figure 1. The figure indicates four cells labeled Ki, i = 1,2,3,4. The star 
in the center of cell K 2 indicates that it is of type A. The star in the center 
right of cell ki indicates that it is of type E. The star in the upper right 
corner of cell K 3 indicates that it is of type NE. The multiple stars in K 4 
indicate that it is of type W, A, or E. 


We begin by labeling the cells according to the behavior of the associated K-equation on 
the cell. In particular. 


(18) 


K{i,j) is of type 


N 

if 

< 

4>i(a) 

< f.i+1 

and rjj+i < 4>2 (k) 

NE 

if 


< 

k) and 

r]j+i < 4>2 (k) 

E 

if 


< 

k) and 

r/j < 4'2 (k) < Vj+i 

SE 

if 

?fi+l 

< 

k) and 

^2{k) < Vj 

S 

if 

e* < 

4>i(k) 

< Ci+l 

and 4>2 (k) < Vj 

SW 

if 


and <h 2 (^) < Vj 

w 

if 


and r]j 

< 4>2(k) < Vj+1 

NW 

if 


and r/j+i < <h 2 (K) 

A 

if 

< 

4>i(k) 

< f,i+l 

and Vj < ^ 2 {k) < Vj+i 


We remark that assumption (|13p implies that every cell k is of the type indicated above. 
A cell of type A is called an attracting cell] a cell of type N, E, S, or W is called a focussing 
cell] and a cell of type NE, SE, SW, or NW is called a translating cell. Furthermore, for 
bookkeeping purposes we find it convenient to place a star in the cell to indicate its type 
(see Figured]). 

As indicated in the introduction we use the switching system to construct the state 
transition diagram. The boundaries of the cells are used in the definition of the set of 
vertices of the state transition diagram. The precise definition is as follows. 


Definition 3.2. Let k = K{i,j) € 1C. The set of faces of k is denoted by V(k) and consists 
of 
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and, additionally, 


Wi^ if K is of type A. 


The complete set of faces is given by 


V := U V(n) 

kGK 

We remark that elements of V(k) have two interpretations: faces of k and hence subsets 
of [0,oo)^, and vertices as elements of the state transition diagram. Whether we are 
employing them as topological or combinatorial objects should be clear from the context. 


Definition 3.3. Given a cell k = K{i,j) we label the faces V(k) as follows 




IS an 


entrance face 
absorbing face 


if $i(k) > 
if ^>i(k) < 




IS an 


entrance face 
absorbing face 


if $i(k) < 
if ^>i(k:) > .fj+i 




IS an 


entrance face 
absorbing face 


if ^>2 (n) > rjj 
if ^>2(«:) < r]j 


Vi 




IS an 


entrance face 
absorbing face 


Wk is an absorbing face. 


if ^>2 (n) < rji+i 
if ^>2(«:) > rji+i 


Let Ve{n) and Va{n) denote the entrance and absorbing faces of V(k). Observe that 


Ve(N)nVa(«) =0. 


As the following proposition indicates, given a pair of adjacent cells that share a face v 
there are constraints on the possible cell types. Using the notation of Definition 13.21 we can 
assume v = j = {^j} x {r]j,rjj^i) or v = v^ j = x {rjj}. In an abuse of notation 

we write 


u.jGH- if gH- 

and 

Uij G if r]j G H" 

for n = 1, 2. 


Proposition 3.4. The following figures show the possible types of two cells that share a 
face V where it is indicated whether n G H"' or n G S” for n = 1, 2. 
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At 
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Ar 

(^) 


(a) 


(in) 


(iv) 


(v) 











(vi) (vii) (via) 


Proof. The arguments for all the cases are essentially the same so we only provide explicit 
proofs in two cases. 

(i) Without loss of generality assume the lower cell is K{i,j). Then the adjacent cell is 
K{i,j + 1). The stars indicate that K{i,j) is of type NW, N or NE. This implies that 
<h 2 («:(i,j)) > By assumption v G Thus d> 2 (K(i,j + 1) = <l> 2 (K(i, j)) > Vj+ii s-iid 

therefore, K{i,j + 1) is of type NW, N, NE, W, A, or E, as indicated. 

(vii) Without loss of generality assume the left cell is Then the adjacent cell 

is n{i + l,j). The stars indicate that n{i,j) is of type W, A or E. This implies that 
Vj < ^ 2 {K{i,j)) < By assumption v G Thus <h2(fi;(i + l,j) = <l>2(K(i, j)), and 

therefore, K{i + l,j) is of type W, A or E. □ 



(ix) 



(x) 


4. (5-Constrained Continuous Switching Systems 

As indicated in the introduction the goal of this section as follows; given a fixed switching 
system, S = S(r, A, H^, H^), construct f^^'^ to define an associated (5-constrained 

continuous switching system (jl]). 
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Figure 2. The shaded regions indicate the four tiles G^{i,j), 

G^{i,j), and G^ = G°(vr) where vr = 

Starting with the observation that 

(19) A = A(S) := imini min - ^i_i} , min I , 

2 1^1=1,...,/ J 

is half of the minimal width of any cell in E we choose 

0 < 5 < A. 

For each k = K{i,j) = (Ci, ?i+i) x {r]j,T]j+i) e 1C, i = 1,..., I, j = 1,..., J define 
G‘^{k) = G‘^{i,j) := - <5] X [rjj + 5,r]j+i - 5] 

■= (0,6 - 5] X (0,r/i - 5] 

■= [6 + <5,6+1 -S] X (0, rji - (5] 

:= (0,6 -6) X [rjj + 6^i+i - <5]. 

For i = 1,..., / and j = 1,..., J define 

G^ivij) = G\i,]) := [6 - <5,6 + <5] X [rjj + 6??i+i - <5] 

:= [6 + <5,6+1 - <5] X [rjj - 5, rjj + 5] 

■= [6 - <5,6 + <5] X (0,?7i - 6] 

■= (0,6 - <5] X [r]j - 5, rjj + <5] 


( 20 ) 


( 21 ) 


G2(0,0) 

G\i,0) 


= G^{i,j) 
= G\i,Q) 
G\v-^^j) = G\Q,j) 


and let 

( 22 ) 


:= [6 - 6 6 + <5] X [rjj - 6, rjj + 6] 


where vr G 11 as defined in ([9]) . We refer to these compact sets shown in Figure [2] as tiles, 
and more precisely G*, f = 0,1, 2 is called an i-tile. 
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We define the continuous nonlinearity : (0, oo)^ —)• (0, oo)^ in steps. Observe that A 
is constant on any given 2-tile G^, and hence, A(G^) is a unique well defined vector. Thus 
we define : Uij G'^{i,j) (0,oo)^ by 

(23) := A{G‘^{i,j)) x G j). 

To define the action of on the 1-tiles we consider four cases. 

• If G then Ai(k(z - l,j) / Ai(K(i, j)) and A 2 (K(i - l,j)) = A 2 {K{i,j)). Thus 
we define 

(24) f2\x) ■= ^2{K{i,j)) X G G^{i,j). 

and choose to be a continuous function on G^{i,j) which agrees with as 
defined by (1^ on G‘^{i — 1 , j) 0 G^{i,j) and G^{i,j) 0 G‘^{i,j) with the constraint 
that for X G G^{i,j) 

(25) min{Ai(K(f - l,j)), Ai(K(f, j))} < fi\x) < max{Ai(K(i - 1 , j)), Ai(K(i, j))} . 

• If G then Ai(K(i - l,j) = Ai{K{i,j)) and A 2 {K{i - l,j)) 7 ^ A 2 (K(i,j)). Thus 
we define 

(26) fi\x) := Ai{K{i,j)) x G G^{i,j). 

and choose to be a continuous function on G^{i,j) which agrees with as 
defined by (f^ on G‘^{i — 1 , j) n G^{i,j) and G^{i,j) n G‘^{i,j) with the constraint 
that for x G G^{i,j) 

(27) niin{A 2 (K(f - 1, j)), A 2 (K(f, j))} < / 2 ^^(x) < niax{A 2 (K(f - 1, j)), A 2 (k(z, j))} . 

• If r]j G H^, then Ai{K{i,j — 1) 7 ^ Ai{K{i,j)) and A 2 {K{i,j — 1)) = A 2 {K{i, j))- Thus 
we define 

(28) f2\^) ■= ^2{K{i,j)) X G G^{i,j). 

and choose to a continuous function on G^{i,j) which agrees with as 

defined by (f^ on G‘^{i,j — 1 ) n G^{i,j) and G^{i,j) n G‘^{i,j) with the constraint 
that for X G G^{i,j) 

(29) niin{Ai(K(f, j - 1 )), Ai(K(f, j))} < f[^\x) < niax{Ai(K(f, j - 1 )), Ai(k(z, j))} . 

• If r]j G H^, then Ai{K{i,j — 1) = Ai{K{i,j)) and A 2 {K{i,j — 1)) 7 ^ A 2 {K{i, j))- Thus 
we define 

(30) fi\x) := Ai{K{i,j)) X G G^{i,j). 

and choose to a continuous function on G^{i,j) which agrees with as 

defined by (f^ on G‘^{i,j — 1 ) n G^{i,j) and G^{i,j) n G‘^{i,j) with the constraint 
that for X G G^(i,j) 

(31) niin{A 2 (K(f, j - 1)), A 2 (K(f, j))} < / 2 ^^(x) < max{A 2 (K(f,i - 1)), A 2 (k(z, j))} . 
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Figure 3. The shaded triangular regions indicate the wide chip C^{K.{i,j — 

1 ), Uj tt) and the narrow chip Observe that 

C^{K{i,j),vi j, 7 r) C G^{vij) audits edges lie in G‘^{K{i,j)) and G^(vr) where 
TT = 

At this point we have defined : (0, oo)^ \ (U^Gn“t(G°(7r))) ^ (0 , cx))^. For each 
TT G n we extend to G'^( 7 r) continuously. 

5. Cells, Chips, and Transversality 

As described in the introduction we construct the lattice of trapping regions using tiles 
and chips. Tiles are introduced in the previous section. In this section we define chips, 
which are closed right triangular regions, and prove results concerning the transversality of 
the vector field of a 5-constrained switching system on the edges of tiles and the hypothenuse 
of chips. The transversality results are used in Section [7] to verify the construction of 
trapping regions. 

A chip is a closed right triangular subset of a 1-tile and is uniquely identified by 

the 2-tile and the 0-tile ^^(vr) that it intersects. With this in mind we denote a chip 

by 

C'”(n, u, 7 r) or C'"'(k, u, 7 r) 

depending on whether it is a narrow chip or wide chip, respectively. If u = or v = vij, 
then the lengths of the edges of a narrow chip are 6 and — ??j )/2 or (^j+i — ^i)/ 2 , 

respectively, while the lengths of the edges of a wide chip are 26 and ( 17^+1 — rjj)!2 or 
— ^i)/2, respectively. Representative wide and narrow chips are shown in Figured 

As indicated above, we use the tiles and chips to construct trapping regions for 6- 
constrained switching systems. For this we need to know that the vector field associated 
with the (5-constrained switching system are transverse to the appropriate edges of the tiles 
and chips. We catalogue this information in the following sequence of propositions. 
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For the remainder of this section S = S(r, A, H^, H^) is assumed to be a fixed 

switching system. Choose 5 > 0 such that 

(32) 3 < 3.(A. 7) min I 3 ,) ■ \/W} ’ 

where /r = /r(S), A = A(S), p = / 9 (S), and 7 = 7 (S) are defined as in (fT4|l . (fT^ . (fTHI) . and 
m, respectively. Note that the above condition implies 5 < A, as 7 < 1 and p < p- 

Proposition 5.1. Consider a € /C. If v G Vein), then the vector field (uni) is transverse 
in to G‘^{k) along the associated edge. 

Proof. Let k = K{i,j)- Assume v = v{i,j). Since v E Vein), n is of type N, NE, A, E, S, 
or SE and <hi(K) > 1,1 + P > Thus the vector field (fT 0 |) is transverse in to G‘^{k) along 

the edge + <5} x [pj + 5, Pj+i — 5]. The arguments for the remain faces are similar. □ 

Similar arguments lead to the following propositions. 

Proposition 5.2. Consider G^{vi j). 

(i) If vj j E Ve(R(i,j)), then the vector field (fTOjl is transverse in to G^^v^j) along the 
edge + x {Pj + <5}- 

(a) If v-^ ■ E Ve(K(i — 1 , j)), then the vector field (fT0]l is transverse in to along 

the edge - <5, x {pj + (5}. 

(Hi) If ^p j^l E Ve{K{i, j))j then the vector field (fTOjl is transverse in to along 

the edge + (5] x {pj+i - 5}. 

(iv) If'Vjpri j_^_i E Ve{K{i — then the vector field (fT0]l is transverse in to 

along the edge — 5,^*] x {pj+i — (^}. 

Proposition 5.3. Consider G^{vi j). 

(i) If E Ve{K{i, j)); then the vector field (fTOjl is transverse in to G^(ip j) along the 
edge + 6} x [pj,pj + <5]. 

(a) E Ve(K(i,i — 1 )), then the vector field (fTOjl is transverse in to along 

the edge + (5} x [pj — 5, pj]. 

(in) If Vi^i j E Ve(R(i,j)), then the vector field (fTU|l is transverse in to along 

the edge {^i+i - (5} x [pj,pj + 5]. 

(iv) If E Ve{n{i,j — 1)), then the vector field (fTU|l is transverse in to G^{vi j) 

along the edge {^i+i — 5} x [pj — 5,pj]. 

Proposition 5.4. Consider G^{i,j). 

(i) If K{i,j) is of type NW, W, or SW, then the vector field (fTU|l is transverse in to 
G^{i,j) along the edge + <5} x [pj,Pj + 5]. 

(a) If K{i,j) is of type SW, S, or SE, then the vector field (fTUD is transverse in to 
G^{i,j) along the edge + (5] x {pj + 5}. 

(in) If K{i — l,j) is of type SW, S, or SE, then the vector field (fTOl) is transverse in to 
G^{i,j) along the edge — <5,x {pj + (5}. 
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(iv) If K{i — l,j) is of type NE, E, or SE, then the vector field (fTOl) is transverse in to 
G^{i,j) along the edge — 5} x [rjj,rjj + c^]. 

(v) If K{i — l,j — 1) is of type NE, E, or SE, then the vector field (fTOll is transverse 

in to G^{i,j) along the edge — (5} x [r]j — S,r]j]. 

(vi) If K{i — l,j — 1) is of type NW, N, or NE, then the vector field (fTOll is transverse 

in to G^{i,j) along the edge x {r]j — <5}. 

(vii) If K,{i,j — 1) is of type NW, N, or NE, then the vector field (fTOll is transverse in 
to G^{i,j) along the edge + <5] x {tjj — 5}. 

(via) If K{i,j — 1) is of type NW, W, or SW, then the vector field (fTOl) is transverse in 
to G^{i,j) along the edge + (^} x [rjj — 6,r]j]. 

We now give a proposition which shows the transversality on the hypotenuse of a chip. 
Consider a 5-constrained continuous switching system ([5|) associated with a switching sys¬ 
tem S = S(r,A, and suppose a narrow chip C^{k, v,7r) or a wide chip 

C'^{k, v,7r) is introduced for S. Recall that 5 > 0 is chosen to satisfy the condition (l3^ . 

Without loss of generality, one can assume, by applying rotation if necessary, that a chip 
C appears associated with k = v = vr = Let H be the hypotenuse 

of the chip G. Under these circumstance, the transversality results can be formulated as 
follows; 

Proposition 5.5. (i) If G = G^{K,v,7r) for n = K{i,j), v = ir^j, vr = {f,i,r]j), k = 

K{i,j) is of type A, N, NE, E, and td = K{i,j — 1) is of type NW, then the vector 
field dH) is transverse in to C along the hypotenuse H. 

(a) If G = G^{k,v,-k) for K = K{i,j), v = vr = j) is of type A, 

N, NE, or E, and td = K{i,j — 1) is of type N or NE, then the vector field (jH) is 
transverse in to C along the hypotenuse H. 

To prove the above proposition, recall that the vector field ([1]) is rewritten as 


X = -Tx + f^^\x) = -r(a: - 


X = — 


From the above definition of f^^\x), we see that = $(k) for any x G G‘^{k) with 

K = K{i,j) or K,{i,j — 1). We first show that the transversality of the vector field ([4]) on i/ is 
reduced to the transversality of an affine vector field. Let p{t) = {l — t)Ho + tIIi (0 < t < 1) 
where Hq and Hi are the end points of H, and let 0{t) is the angle (measured counter¬ 
clockwise from the positive direction of the xi-axis) of the vector V{x) = —Fx -|- f^^\x) 
evaluated at p{t). In other words, 



Observe that the types of k, td and the definition of f^^\x) on G^{v) imply that the angle 
9{t) must satisfy 0 < 9{t) < vr for any t G [0,1]. 
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V{p{t)) ■ V at any point of p(t) G H {t G [0,1] is positive. Let 6*min be the minimum of 6{t), 
and let tmln £ [0, 1] be such that 0(tmin) = dmm- 

Lemma 5.6. Let V^am = b^(p(imin))- If^mm • u > 0, then V{p{t)) • u > 0 for any t G [0,1]. 


= R{t)r cos{0{t) — (p), 


Proof. Observe that 

where —tt < O^m — P < d(t) — p < 7r/2. If 0 < 6(t) — p < 7r/2, then we immediately obtain 
the positivity of V{p{t)) ■ u. So we consider the case —tt < 0min — P < d(t) — p < 0. In this 
range, the cosine function is monotone increasing, and hence cos(0min — p) ^ cos{6(t) — p). 
Therefore, if I4iin • u > 0, then we have 0 < cos(0min — p) < cos{6{t) — p), and therefore 
V{p{t)) • u > 0. □ 

Now define <I>* = —r“^I4iin and consider the affine vector field V*{x) = —r(x — d)*). 


Observe that, from the definition of on G^(u), <I>* = 


satisfies 


for any ^ G'B, 77 G H. Similarly, it also satisfies 

p<M-v\ 

for rj = rjj. 

In order to prove Lmin • u > 0, we first consider the narrow chip C = C^{k, v, tt). 
In this case, we can use coordinate system with origin at {f.i,Pj) to evaluate 


pit) = (1 -t) +t 


a/2\ _ f6 + (a/2 — 5)t 


and one can choose u = (^^^2 ^ ’ ^bere a = — ^i(> 2A). Therefore the inner product 

V*{p{t)) ■ V defines a function Tit, 6) which is affine in t G [0,1] and quadratic in 5. More 
explicitly, we can obtain Tit, 6) = Kit)6‘^ + L(t)(5 + M(t), where 


(33) 

Kit) = 

(1 - t)7i + t72 


(34) 

Lit) = 

(71 -72)(a/2)t-7i(^i -C 

1 

to * 

1 

(35) 

Mit) = 

( 0 / 2)72 (4>2 -7j) 



In fact, M does not depend on t, M(t) = M, and is strictly positive. Also Kit) is strictly 
positive because 71,72 > 0 and t G [0,1]. So the quadratic function Tit, 5) = Kit)5‘^ + 
Lit)5 + M is always positive for <5 > 0 if the discriminant D = Dit) = L(t)^ — 4iL(t)M 
is negative, or if Z) > 0 and L > 0. So the only case where Tit, 6) can take a negative 
value is the case D > 0 and L < 0. In this case, we can easily see that T(t, d) > 0 for any 
S G (0,(5_(t)), where 5-(t) is a smaller root of the quadratic function Tit,5) = 0, namely 


5. it) = 


-Ljt) - 
2Kif) 
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In order to obtain an estimate which does not depend on t, we observe that 6-{t) > 
—M/L{t) in case L{t) < 0. This easily follows by comparing the graph of T{t,6) and its 
Substituting above and using the estimates (fT^ . (fTSll . (fT6|) . (fT9|) . we 

_(a/2)72(^2 -^i)_ 

-(71 - l 2 ){a/ 2 )t + 7i($* - + 72(d'^ - rjj) 

> _ (a/2)72(^2 - Vj) _ 

l7i -72 |(o/ 2) + 71 -^i) + 72(^2 -Vj) 

_1^2-^jl_ 

1 ( 71 / 72 ) - 1| + (2/a){(7i/72)|$t - + \^*2 - r]j\} 

> _ ^ _ 

- 1(1/7)-l| + (l/A){(l/7)p + p} 

_A//7_ 

(1-7)A + (1 + 7)/9 
^ A/i7 
A + 2/3 

Therefore, in the case of the narrow chip, if we choose 6 so that 0 < <5 < + 2p), 

then we can conclude that the vector field V*{x), and hence V^^\x) as well, is transverse 
in to the chip along its hypotenuse. 

For the wide chip, we can argue in the same manner, but the estimates become more 
complicated. Using 


and 12 = 

( ) 

\a/2-6j 

, one can similarly define 



r(t, ,5) := V*{p{t)) • p = + L{t)6 + M 

where 



(36) 

K{t) 

= 271 - 72 - 2(71 - 72)t, 

(37) 

L{t) 

= (71 - 'l 2 )at + (0/2)72 - 271(^1 - ^i) - 72(^2 

(38) 

M{t) 

= (a/2)72($2-?/j)- 


As in the case of narrow chip, M does not depend on t and is positive. However, K{t) can 
change its sign. In case K(t) > 0, the same argument works, and we obtain that —M/L{t) 
gives a bound for 6 in case L{t) < 0, otherwise T{t, <5) > 0 for any d > 0. 

Observe that K{t) vanishes at some to £ [0,1] only when 72 > 271. In case K{to) = 0, 
T(to,<5) = 0 when 5 = —M/L{to) and 

L{to) = (7i-72)aio+(a/2)72 - 27i(^t-^i)-72(^2-^i) = (7i-72)a-27i(^>i-^0-72(^2-^i) 


tangent line at 5 = 0. 
obtain 

M 

~W) 









2IDMAS GEDEON, SHAUN HARKER, HIROSHI KOKUBU, KONSTANTIN MISCHAIKOW, AND HIROE OKA 


SO the upper bound of 6 for T{tQ,5) > 0 is <5 < —M/L{to) which is the same as before, or 
more precisely, 

M ^ _ { a / 2 ) j 2{<^*2 - Vj ) _ 

L{to) (72 - 7i)a + 27 i(^>| - Ci) + 72(^2 “ 

_1^2 -^jl_ 

2|1 - (71/72)1 + (2/0)12(71/72)I ^>1 - Ci \ + \^2 - Vj \} 

_A/17_ 

2(1-7)A + (2 + 7)p 
^ A/i7 
“ 2X + 3p' 

In case K{t) < 0, which occurs when 72 > 271 and 0 < t < to; we have T{t,5) > 0 if 
0 < 5 < (5+(t) = {—L{t) — yjD{t))/2K{t), where S+{t) is a larger root of the quadratic 
equation T{t,6) = 0. Using an inequality 

y/x + h > "s/x + ^ ^ : (x > 0, /i > 0) 

2 vx + h 

we have, independent of whether L{t) > 0 or not, 


5+{t) > 


L { t ) + \ L { t )\ + 


1 4(-K(I))M 

2 y/Lit)^+4{-K{t))M 


2{-K{t)) 


> 


M 

VU(t)2 + 4(-if(t))M’ 


Since L{t)‘^ + A{—K{t))M < max{2L(t)2,2 x A{—K{t))M}, we finally have an estimate 


<5+(t) > min 


M 1 I M \ 
V 2 \ L { t )\' - K { t)j 


The first term in the min can be treated exactly the same way, up to the constant l/\/2- 
The second term can be treated as follows: 

M ^ (0/2)72(^^ -r/j) ^_ (o/2)(^* -rjj) _ ^ A/17 ^ 

- K { t ) 1271 - 72I + 2|7i - 72I 12(71/72) - 1 | + 21(71/72) - 1 | (2 -7)+ 2(1 -7) 

Putting all the above estimates together, for both the narrow chip case and the wide 
chip case, we define 

(39) \/M|. 

The above argument proves that, for any 6 with 0 < <5 < d*, the vector field U*(x), 
and hence U^'^^(x) as well, is transverse in to the (both narrow and wide) chip along its 
hypotenuse. This completes the proof of Proposition 15.51 




























GLOBAL DYNAMICS FOR STEEP SIGMOIDAL NONLINEARITIES IN TWO DIMENSIONS 


21 



_ 


%i±l 


_ 


D,J( 


Wk 

'^i,j c 



) D+l,j '^i,j ( 













Fj 



(a) (b) (c) 

Figure 4. The directed graphs or multivalued maps Ve(K) F Vain) 
for K = (a) K is attracting and hence of type A. (b) k is focusing of 

type N. (c) K is translating cell of type NE. 

6. State Transition Diagram 

Definition 6.1. Given a cell k the directed n-graph, F^- Vei^) ^ Vain), is dehned by 

V G F{u) if and only if v G Va{i^) and u G Vei^). 

In addition, if G V(n), then F^iw^) = w^- 

To give some geometric perspective to these definitions the reader is referred to Figured] 
where the directed graphs F^ are shown for various types of cells k. The following propo¬ 
sition, whose proof is left to the reader, relates the definitions of entrance and absorbing 
faces to the dynamics of the associate K-equation. 

Proposition 6.2. Let v G Ve{i^)- For every x G v, there exists a unique t = t{x) > 0 and 
a v' G Vain) such that 

x) G cl{v') or lim ijjKit,v) = ^{k). 

t^OO 

In the latter case, k is an attracting cell and Fi^{v) = w. 

Definition 6.3. Given a switching system S the associated state transition diagram 
T”: V ^ V is the directed graph with vertices 



where ~ is equivalence relation v ^ w when v G V(«:), w G V{k') and kH k' intersect in a 
1 -dimensional interval that correspond to both v and w. 

We label the possible configurations of and F^' under the assumption that k and k' 
share a common face. 
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Definition 6.4. Consider the state transition diagram T:V^V for a switching system. 
Let i; E V. li v = Wij E V{K{i,j)), then v is called an minimal vertex. Assume v E 
V(k) n V(k'). If there exists u E V{k) and w E V{k!) such that u ^ v ^ w, then u is a 
transparent vertex. If there exists u E V(k) and w E V(k') such that u — >■ u and w ^ v, 
then u is a black vertex. If there exists u E V(k) and w E V(k') such that u —>■ u and v ^ w, 
then u is a white vertex. The sets of minimal, transparent, black, and white vertices are 
denoted by A4, T, 13, and W, respectively. 

Proposition 6.5. Consider a state transition diagram T'.V^V for a switching system. 
The minimal, transparent, black, and white vertices partition V. 

Proof. We need to demonstrate that AiCTcBCW = V and that M., T, B, and W are 
mutually disjoint. By definition Ad n (TU B U W) = 0. Thus we can restrict our attention 
to vertices that are associated with faces of cells. 

Observe that from Figure|4]if v E V(k), then there is either an edge to v or an edge from 
V in JCc; but not both. The first statement follows from the existence of the edges. The 
second statement follows from the fact that it is not possible to have both types of edges 
within one cell. □ 


For the remainder of the paper we make the following assumption 

P: The state transition diagram T :V does not contain a black vertex. 

We hasten to add that this is not an unreasonable assumption. To understand why, 
consider the following result that is easily checked. 

Lemma 6.6. If v{i,j) is a black vertex, then 

Ai{K{i - l,j)) > 'jifi > Ai{K{i,j). 


Observe that starting with a regulatory network this implies that xi has a self-edge that 
corresponds to repression. For many biological applications this type of self-regulation is 
better modeled by two nodes xi and y where xi activates y and y represses xi. 

From a mathematical perspective. Lemma 16.61 can be used to show that given a system 
of the form m it is possible to approximate the nonlinear functions /i using piecewise 
constant functions in such a way that the state transition graph f or resulting switching 


system does not co ntain black walls. This approach is discussed in lEdwards et al.l [201 


and in future work ICummins et al.l 2015al | . 


We conclude this section with three results concerning the structure of forward invariant 
sets under T. 


Proposition 6.7. Consider M E lnvset^(V, T") under assumption If v G Af^, then 
V G Ai UT. Furthermore, if v gT, then there exist distinct u G M and w G such that 
u ^ V ^ w. 


Proof. By Proposition 16.51 and assumption P, it is sufficient to show that v 0 W. If u E W, 
then F~^{v) = 0, hence v ^ a contradiction. 

Since v G by definition there exists u G F~^{v) CAf. If u E T, then v does not 
have a self-edge hence u ^ v. By definition, if u E T, then there exists w G V such that 
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Figure 5. Possible intersections of with a translation cell. Vertices of 
are indicated as solid circles. 

V ^ w. This implies that w G J^(AA). By definition of a forward invariant set C M 

and hence w G J\f. Since v ^ w, w £ □ 

Proposition 6.8. Let Af G Invset^(J^). IfV{K) CiAf^ n T 7 ^ 0, then Ve(N) H AA / 0 and 
Vaitd) n AA c AA°. 

Proof. Let v G V(n) ClAf^ Pi T. By Proposition 16.71 there exist exist distinct u £ Af and 
w £ Af^ such that u ^ v ^ w. li v £ Ve{tf), then Vai^) = and hence Vai^) C Af^. 

li V £ Va(fi;), then u £ Af £\ Ve(N) and J^niu) = Va{n). Therefore, Vain) C Af^. □ 

Proposition 6.9. Consider Af £ lnvset'’'(J^) under assumption If v £Af^nV{K), then 

Vain) C AA°. 

Furthermore, if Af^ n V(n) 7 ^ {wk}, then Vein) n A^ 7 ^ 0. 

Proof. First, assume v £ Ve{iC). Then, Va(N) C Tk{v), and hence, by Proposition 12.51 

Va(K) C A. 

Now assume v £ Vain). To indicate the line of reasoning we provide the proofs of two of 
the three cases. Assume k is an attracting cell. This implies that Va{if) = {r^k} = {u} C 
Va{iCj. Assume k is a translating cell. By Proposition 16.71 there exists u £ Af and w £ Af^ 
such that u ^ V ^ w. Since v £ Va{K), by P rt G V(k). Furthermore, u ^ v implies that 
u £ Vein). Thus by Proposition 12.51 .F.- (u) = Va{K) C Af^. □ 

7. Constructing Trapping Regions for (5-Continuous Switching Systems 

Let Af £ Invset’^(J^). Our goal is to identify a compact set Nj^ C (0, 00 )^ such that Nj\f 
is a trapping region for 

(40) X = —Px + A^‘^^(a:). 

The region Nj./ is constructed using the above defined tiles and chips. Furthermore, as is 
made clear below the identification of Nj\f is done locally. 
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Definition 7.1. For each {i, j), i = 0,..., I, j = 0,..., J the associated elementary domain 
is defined to be 

E{i, j) ■■= {Ci-i + S, Ci+i - (5) X {rij_i + 6, rjj+i - 5) . 

The set of associated vertices is denoted by £{i,j) C V and is defined to be the union 
of vertices for which the associated face v satisfies v (1 E{i,j) / 0 or a face where 
K n E{i,j) / 0. We use £{i,j) to denote the union of vertices in £{i,j) and the vertices 
associated with the faces of the cells in x [r/j-i, r/j+i]. 

We leave it to the reader to check that 
(41) (<5,00)2= y 

z = 0, 

j = 0,...,J 

Fix an elementary domain E{i,j). In the applying the following rules, unless otherwise 
specified, it is assumed that the cells k, the faces v, and vr = rjj) G 11 all intersect E{i,j) 
nontrivially. 

Rule 0: If Wk G n A4 n £{i,j), then 

G2(k) c IYaT- 

Rule 1: Let u G AA^nTnT(i, j) and let r, s G £{i,j) such that r ^ u ^ s. Consider 
the grid elements Kr and Ks such that r, n G V{nr) and u,s G V(ks). Then 

G'^{Kr) U G‘^{ks) U G^(u) C Nj^. 

Rule 2: Assume 

G‘^{tia) U G‘^{Kj3) U G‘^{Kj) U G^{Val3) U G^{Vj3^) C N_\f, 
where Vaj 3 = Kq n and = upCi k^. If Vay G Vai^y) or G Va(N/ 3 ), then 

G°(7r) C Nj^. 

Rule 3: Assume G^(k) U G^(uq) U G°(7r) C AA/" where Uq G V(k) n £l(f,j). IfV(K)n 
£ii,j) C Vein), then 

G'^iK,V0,Tr) C AAt 

where U /3 G V(k) n T(i, j). 

Rule 4: Assume G^in) C Nj^ and G'^ina) Nf,/. Let = V{Ka) nT(f,j). 

If Vai3 0 Vain), then 

C inj, Vci'y,'^) <G Elj\f 

where = Kq n K-y. 

Rule 5: Assume G'^ina) U G^(k^) U G°(7r) C Nj\f and Vap = Ka H np. If Vina) H 
£ii,j) C Veina) and Vinp) n £ii, j) C Veinp), then 

G\vap)^Nj,f. 
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Definition 7.2. For each Af E Invset'^(J^) define Nj^ C (0, oo)^ to be the union of the 
minimal collection of tiles and chips that satisfies Rules 0-5 over all elementary domains 
i = 0,...I, j = 0,...,J. 

For future reference we highlight the following remark. 

Lemma 7.3. If Af = (/i, then Nj^ = A ^0 = 0. 

Proposition 7.4. If Af ^ Invset'^(J^), then Nj^ E AN bhd ((/>). 

The goal for the remainder of this section is to prove Proposition 17.41 We will do this 
by proving that for each Af E Att(J-'), Nj\/ is a trapping region. This in turn is done by 
considering all possible forms of intersection of Nj^ with all possible elementary domains 
E(i,j) and checking for transversality as we proceed. For this we make use of the following 
objects. 

Definition 7.5. Let e denote a boundary edge of a 2-tile G^(k), a 1-tile G^(n), or a 0-tile 
G^{i,j). We say that e is interior to E{i,j) if 

enE{i,j) / 0 

and exterior to E{i,j) if e is not interior to E{i,j), but 

e C d{E{i,j)). 

A boundary edge e of a chip is interior to E{i,j) if e n j) 

As the following lemma indicates it is sufficient show that for every E{i,j) along any 
boundary edge of N that is interior to E{i,j) the vector field of (14011 is transverse in with 
respect to N. 

Lemma 7.6. If N is a union of tiles and chips and for every elementary domain E{i,j), 
z = 1,..., j = 1,..., J, the interior edges of N with respect to E{i,j) are transverse in, 
then N is a trapping region. 

Proof. By (j4ip the collection of elementary domains covers the phase space (0, oo)^. By 
definition 6 < X which is half the minimal cell width (see m and thus every point x E 
(0,oo)^ lies in the interior of some E{i,j). Thus if the interior edges of N with respect 
to E{i,j) are transverse in for all elementary domains, then every boundary point of N is 
transverse in. Therefore A^ is a trapping region. □ 

To simplify the application of the rules we use the symmetry of the elementary domain. 
By Proposition 16.71 if AA E Invset(T'), then Af^ C A4 U T. Assume that there exists 
V E Af^ n T. This implies that v represents the intersection of two cells ki and K 2 . Since 
n E T, there exists u,w G V such that u ^ v ^ w. Without loss of generality we assume 
that u E V(n 2 ) and w E V(zti). Symmetry, i.e. a reflection or rotation, allows us without 
loss of generality to assume that K 2 = n{i,j) and ki = K{i + 1 , j) and thus we will refer to 
n as a transparent face moving east. We will study the implications of Rules 0 - 5 on the 
elementary domain E{i,j) shown in Figure [ 6 ] where the neighboring cells of interest are 
K 3 = K{i,j — 1 ) and A 4 = K{i + l,j — 1 ). 
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We begin by making simple observations concerning Rules 0-5. If AA^nT r\E{i, J) / 0, 
then Rule 1 implies that 

(42) G^(ki)UG^(k2)UG\v) cN^f 
as indicated in Figure [6Kb). 

Rule 2 determines if G^ tiles belong to Nj\f. Therefore, a necessary condition for 
G^(i,j) C Nj\f is that at least three of the four G^ and at least two of the four G^ tiles in 
the elementary domain E{i,j) belong to N_\f. 

Observe that Rule 4 is only applicable if given an elementary domain E exactly three 
of the four associated G^ tiles belong to N_^f. Furthermore, if Rule 3 implies that 
C^{k,Vi 3 , {i,j)) C Nj^ and Rule 4 implies that G'^{K,vp, {i,j)) C N_\f, then we can ignore 
Rule 3 since 

G^{K,Vf^,{i,j)) C G^{K,Vf^,{i,j)) C N^f. 

Observe that if Rule 5 applies then, Rule 3 applies to ki and K4. However 

G"-{Kn,vi,Tr) C G^{vi), n = l,4 
and hence we can ignore Rule 3. 

Proposition 7.7. Let {va^} = V(kq.) O V(k^) where Kn O E{i,j) n = a, ( 3 . If 

G‘^{Ka) <t- then G^{va^) (/L N^f. 

Proof. Only Rule 1 and Rule 5 require that G^{va0) C Nj^. Both these rules are based 
on G‘^{Ka) C Nj^. Thus the minimality of Nj^ implies that G^{va/ 3 ) Nj^. □ 

Proposition 7.8. Let Mo,Ml G Invset^(T'). If Mq =Mi, then 

hlMo = ■ 

Proof. The only rules that directly depend on the elements of M are Rule 0 and Rule 1. 
These rules are given in terms of M^. □ 

Proposition 7.9. If M = {Kk | A: = 1,..., K} C Ml, then M = M^ G Att(T'), 

Nm= U 

k=l,...,K 

and Nj^ is a trapping region. 

Proof. By Rule 0 

IJ G^{Kk) C N^. 

k=l,...,K 

We now show that this is the minimal collection of tiles and chips that satisfy Rules 0-5. 
Given that M Cl T = 0 , the only way to require the existence of a G^ tile is through Rule 
5. However, Rule 5 requires the existence of G^{i,j). The requirement for the existence 
of G^{i,j) follows from Rule 2, which in turn requires the existence of a G^ tile. Thus, 

(43) Nm= U G\Kk). 

k=l,...,K 


GLOBAL DYNAMICS FOR STEEP SIGMOIDAL NONLINEARITIES IN TWO DIMENSIONS 


27 


'k 

lAr-i 

K2 ^ 

V2 

k k 

^ 1 * * 

'k k: k: 

k k lA: ( 

k k k 

A 

'vt * * 

k k k 


(a) 


G 

G‘^{k2) 


v) 

; G2(ki) 

- 


-^ 


•L 

; 


(b) 


Figure 6. Black filled dots indicate vertex that belongs to J \ f ^. Blue filled 
dot indicates that vertex may or may not belong to (a) The vertex 
u is a transparent face moving east. The associated four cells are Ki, i = 

1,..., 4. The possible cell types of ki and K 2 are indicated by the stars, e.g. 
center is A, upper right corner is NE, middle rights is E. Without further 
assumptions there are no restrictions on the cell types of and K 4 . (b) 

Shaded region indicates tiles that belong to Nj^ under the assumption that 

WO n r / 0 . 

For each k = 1,..., W, each face is identified with a vertex that belongs to Ve{nk)- 
Thus, by Proposition 15.11 G‘^{Kk) is a trapping region and hence Nj^ is a trapping region. 

□ 

The following result is an immediate application of Propositions 17.81 and 17.91 
Corollary 7.10. Let J\f G Invsef*'(. //W^ = {na, | A: = 1,... ,K} C M, then 

Nm= U 

k=l,...,K 

and Nj\f is a trapping region. 

For the remainder of the argument we assume that ^ A4 or equivalently that W® n 
T 7 ^ 0. Furthermore, taking advantage of the above described symmetry we always assume 
that V = Ujj G n T is transparent east. 

Proposition 7.11. Under the assumption that v = v^ j G H T the cells that intersect 
E{i,j) must be of the types indicated in Figure\B(a) and as indicated in Figure\^h) it must 
he the case that 

G\ki)UG\k2)UG\v) CNm. 

Proof. The assumption that v = is transparent east implies that <1 >i(ki) > and 
‘hi(At 2 ) > Thus, Ki is of type N, A, S, NE, E, or SE and K 2 is of type NE, E, or SE. 
Rule 1 implies that G‘^{ki) U G‘^{k 2 ) U G^{v) C N_f^. □ 
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Figure 7. Unfilled dot implies that vertex does not belong to (a) 
Possible cells types in the setting of Case 1. (b) Minimal set of tiles in Nj\f 
in Case 1. (c) Possible cell types in the setting of Case 1(a). 
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Figure 8. (a) Cells in Nj^ associated with under the assumptions 

of Case l(a)(i). (b) Possible cell types in the setting of Case l(a)(ii). (c) 
Tiles for Nj^ associated with E{i,j) in Case l(a)(ii). 


Proof of Proposition (1) Assume ui € and vi is transparent north. By Rule 

1, G^(k 4 ) U G'^(ui) C Nj\j-. The possible cell types are indicated in Figure [T]) a). 

(a) Assume V 2 0 By Proposition 16.91 ko cannot be of type SE and thus by 
Figure [7Ka) must be of type E or NE. By Proposition 13 . 4l fil lipj G H^, then ^2 
being of type E or NE implies that K 3 is of type NW, N or NE. Similarly, by 
Proposition I3.4r vl if rjj G H^, then K 3 is of type NE, E, or SE. The possible 
cell types are indicated in Figure [7Kc). Since v 0 Ve(K 2 ), Rule 3 does not 
apply to K 2 . 

(i) Assume V3 0 J\f^ Because of the cell type of ki, Rule 2 does not apply. 
Because G^{Tr) ^ Nj^, Rule 3, Rule 4 and Rule 5 do not apply. We 
claim that we are in the setting of Figure [81(a). In particular, we need to 
argue that G^(k 3 ) ^ N_\f. Observe that only Rule 0 and Rule 1 require 
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Figure 9. (a) Possible cell types in the setting of Case l(a)(iii). Observe 
that the possible cell types of K 2 and K 3 cannot occur simultaneously. Thus 
this case cannot occur, (b) Possible cell types in the setting of Case 1(b). 
(c) Minimal set of tiles for Case 1(b). 


the introduction of a 2-tile. By Figure [71(c) K 3 is not of type A and 
hence Rule 0 does not apply. Assume Rule 1 forces the introduction 
of Then there exists u E such that u E V(n 3 ). Again, by 

Figure [ 7 j(b)u 2 or V 3 belong to VaiKs) and hence V 2 or vs belong to a 
contradiction. Thus, we are in the setting of Figure[ 8 ](a). Proposition l5.ll 
guarantees the desired transversality on the edges of G‘^{k2) and G^(n 4 ). 
Furthermore, Proposition 15.2f i - ii) and Proposition 15.30 - ii) guarantees 
the desired transversality on the edges of G*^( 7 r). 

(ii) Assume V 3 E and V 3 is transparent east. The assumption that V 3 is 
transparent east implies that N4 is of type N or NE and K3 is of type 
NE, E, or SE. By Proposition 16.91 if K3 is of type NE or K2 is of type 
SE, V2 E Af^, a contradiction. Thus, the possible cell types are as shown 
in Figure EKb). 

By Rule 1 G^(k 3 )uG^(u 3 ) C Nj\f. Applying Rule 2 to G^(k 3 )uG^(n 4 )U 
G^(ni) U G^{v3) U G^(ui) implies that G^{Tr) C A(yy. This implies that 
U 3 ^ Ve(N 3 ) and hence. Rule 3 does not apply to K 3 . Rule 4 and Rule 
5 do not apply and hence we are in the setting of Figure (SKc). Propo¬ 
sition IQ guarantees the desired transversality on the edges of G'^{k 2 ) 
and G‘^{ks). Proposition I5.4l iv-vi guarantees the desired transversality 
on the edges of G^{k2). 

(iii) Assume U 3 E Af^ and U 3 is transparent west. Then N 3 is of type NW or 
N. See Figure [91(a). Since U 3 E Af^ Proposition 16.91 implies that V 2 eA^o 
a contradiction. Thus this case cannot occur. 

(b) Assume V2 E Af^ and V2 is transparent south. This implies that K2 is of type SE 
and K 3 is of type E or SE. See Figure[9l(b). By Rule 1, G^(k3)UG^(u2) C Nj\f. 
By Rule 2, G^{tt) C Nj./. See Figure El(c). 
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Figure 10. (a) Possible cell types in the setting of Case 1(b). (b) Possible 
cell types in the setting of Case l(c)(i). (c) Tiles in Nj^ that are related to 
E{i,j) in Case l(c)(i). The vector held is transverse in along the interior 
edges indicated in red. 


(i) Assume 0 Since V 2 G if K 3 is of type E or SE, then 

^^3 G Va(K 3 ), a contradiction. Thus, K 3 is of type W, A, SW or S, and 
hence, by Proposition I.S.dI x) K 4 is not of type N or NE. Therefore, 
the set of possible types is as in Figure fTOl ai from which we note that 
V'i G Va(K 4 ). Since vi G by ProDosition l6.9l uq G Af^, a contradiction. 
Therefore, this case cannot occur. 

(ii) Assume V3 G Af^. By Rule 1, G^(u 3 ) C Nj^. Combining this with 
the information from Figure [9Kc) we observe that there are no interior 
boundary edges to check. 

(c) Assume V 2 G Af^ and V 2 is transparent north. 

(i) Assume V 3 0 Af^. By Proposition 16.91 the assumption that V 3 ^ Af^ 

implies that K 3 is of type NW or N and K 4 is of type N or NE (see 
Figure [TOl bB. By Rule 1 U G^{v2) C Nf^. By Rule 2 G°( 7 r) C 

Nj^f. Rule 3-5 do not apply, thus we are in the setting of Figure fTOl cb 
Thus, by Proposition lS. li the vector held is transverse in along the bound¬ 
ary edges of G^(a 3 ) and G^(a 4 ). By Proposition I5.4lf vi-viii the vector 
held is transverse in along the boundary edges of G*^( 7 r). 

(ii) Assume vs G Af^. By Rule 1, G^{v 3 ) C N_\f. Observe that there are no 
interior boundary edges to check. 

(2) Assume vi G Af^ and vi is transparent south. Recall that Figure [61(a) provides an 
upper bound of the types of cells. Since vi is transparent down, vi G Va(Ki) and 
vi G Ve{Ki). Therefore, ki is of type S or SE and K 4 is of type IT, A, E, 51T, S 
or SE. The possible types of K 4 preclude the possibility of K 3 being of type NW 
or N. Thus we are in the setting of Figure fTR a). 

By Rule 1, G‘^{k4)uG^{vi) C Nj\f. By Rule 2, G^{7r) C Nj^. See Figure fTlT b). 
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Figure 11. (a) Possible cell types for Case 2. (b) Minimal set of grid 
elements in Nj\f related to E{i,j) in Case 2(a). (c) Possible cell types for 
Case 2(a). 


(a) Assume V2 0 Figure fTlT al provides an upper bound on the cell types. 
By Proposition 16.91 ko is of type E or NE. This in turn, by Proposition 13.41 
implies that must be of type E, NE or SE (see Figure [TlT cB. 

Therefore, v E Va(N 2 ) and hence Rule 3 does not apply to K 2 - Thus we are 
still in the setting of Figure fTlT b). 

(i) Assume V 3 0 By Proposition 16.91 there cannot be an edge —>■ U 3 . 

Thus, A 4 is of type A, S, E, or SE (see Figure [T2]( a)). We claim that 
G^(k 3 ) ^ N_f^. Observe that only Rule 0 and Rule 1 require the 
introduction of a 2-tile. K 3 is not of type A and hence Rule 0 does not 
apply. Assume Rule 1 forces the introduction of G'^{k 3 ). Then there 
exists u E AA^nT such that u E V(k3). Since K3 is of type E, NE or SE, 
V3 E Va(K3) and hence V3 E a contradiction. By Proposition 17.71 
G^{v 2 ) and G^(u 3 ) are not in N_^f. Furthermore, Rule 5 does not apply. 
Thus we remain in the setting of Figure [TTT bh 

(A) Assume K 3 is of type NE. Then Rule 4 does not apply to K 3 . 
By Rule 3 applied to K 4 , C'”(k 4 , U 3 , tt) C N_\f. See Figure [T^ bi. 
We need to check that the interior edges of dNj\f are transverse 
in. Observe that V2 E Ve(K2) and V3 E Ve(K4), hence by Proposi¬ 
tion [Q we have the desire transversality along the edge of G'^{k 2 ) 
and G^(a 4 ). By Proposition l5.11 iv-vii we have the desired transver¬ 
sality for G^ (tt) . By Proposition 15.51 we have the desired transver¬ 
sality for C'"’(a4, U 3 , tt). 

(B) Assume K 3 is of type E or SE. Rule 4 applied to K 3 implies that 
C'“’(k 4 , ^ 3 , tt) C Nj^. See Figure fT^ ci. The desired transversality 
follows from the argument used in the previous case. 
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Figure 12. (a) Possible cell types for Case 2(a) (i). (b) Set of grid elements 
and chips in Nj\f associated with E{i,j) in Case 2(a)(i)(A). (c) Set of grid 
elements and chips in Nj^ associated with E{i,j) in Case 2(a)(i)(B). 
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Figure 13. (a) Possible cell types for Case 2(b). (b) Necessary tiles in Nj^/ 
for Case 2(b)(i). (c) Possible cell types for Case 2(c). 


(ii) Assume V 3 G and is transparent east. If we take this configuration 
and rotate it counterclockwise by 90° then we are in case 1 (c) (i) for which 
we have already shown the desired transversality. 

(iii) Assume V3 G and V3 is transparent west. Figure m^) provides an 
upper bound on the cell types. However, U 3 is transparent west implies 
that AC 3 cannot be of type NE, E, or SE. Thus this case cannot occur. 

(b) Assume V 2 G and is transparent south. Figure [UK a) provides an upper 
bound on the cell types. Since V 2 is transparent south, k ,2 is of type SE and 
K 3 is of type W, A, E, SW, S, or SE as indicated in Figure fTHlf a). 

By Rule 1, G‘^{k 3 )UG^{v 2 ) C Nj^. By Rule 2, G°( 7 r) C Nj\f. Thus Nj\f must 
contain at least the tiles indicated in Figure [THIbb 

(i) Assume V 3 0 Af^. This implies that K 3 cannot be of type E or SE and K 4 
cannot be of type W or SW (see Figure [TSl cB. Therefore V 3 G W. Thus 
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Figure 14. (a) Possible cell types for Case 2(c). (b) Possible cell types for 
Case 3. (c) Possible cell types for Case 3(a). 


Rule 5 applies and G^(i; 3 ) C Nj\f. Observe that there are no internal 
edges against which to check internal tangencies. 

(ii) Assume G By Rule 1, G^{v 3 ) C Nj^. Observe that there are no 
interior boundary edges to check. 

(c) Assume V 2 G and V 2 is transparent north. Figure [TlT al provides an upper 
bound on the cell types. Since V 2 is transparent north, ^3 must be of type 
NE and K 2 cannot be of type SE. The adjacency of K 3 and ^4 implies that 
K 4 cannot be of type W or SW . See Figure [UK a). 

By Rule 1, UG^(u 2 ) C Nj\f. By Rule 2, G^{'k) C Nj^f. The necessary 

tiles in Nj^/ are indicated in Figure [TSKb). 

(i) Assume U3 0 Because V2 G and K3 is of type NE, by Proposi¬ 
tion [619] U3 G a contradiction. Thus this case cannot occur. 

(ii) Assume U3 G A. By Rule 1, G^(u3) C Nj\/. Observe that there are no 
interior boundary edges to check. 

(3) Assume vi 0 Af^. Figure [TKa) provides an upper bound on the cell types. Since 
vi ^ Af^ there is no edge u —>■ ui. Therefore the cell ki is of type A, NE, N or E 
and vi G Ve(Ni). Adjacency of ki and K4 implies that K4 is not of type W or SW. 
See Figure [Tilf bl. 

(a) Assume V 2 0 Af^. Figure [T41f bl provides an upper bound on the cell types. By 
Proposition 16.91 we conclude that K 2 is of type NE or E. By Proposition 13.41 il 
or (iv) the adjacency of K 2 and K 3 implies that K 3 must be of type NW, N, 
NE, E or SE. See Figure [T41f cl. The necessary set of tiles in Nj\j- is given by 
Figure EKb). 

(i) Assume V 3 0 Af^. 

(A) Assume G‘^{ks) Nj^ and G^(k 4 ) N^j-. By Proposition 17.71 

G^{vn) N_^f, n = 1,2,3. Rule 2 is not applicable thus, G‘^(vr) ^ 
Nj^. Thus the set of tiles is not changed from that of Figure [ 6 Kb). 
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Figure 15. (a) Necessary grid elements for Case 3(a)(i)(C). (b) Tiles in 
Case 3(a) (i)(C)(1). (c) Tiles in Case 3(a)(i)(C)(II). 


By Proposition 15.11 the lower face of Gl and by Proposition 15.21 the 
lower left face of is transverse in. 

(B) Assume G^(k 3 ) C Nj^. By Figure fTllc) K 3 is not of type A. Thus 
there exists u E V(k 3 ) H ^ n T. By Proposition 16.81 there exists 
u E Veins) r\Af^. By assumption u E V(k 3 ) \ {v 2 ,vs}- We leave 
it to the reader to check that the given the possible types of K 3 by 
Proposition 16.91 either V 2 or vs is in A, a contradiction. Thus, it is 
not possible for G'^ins) C Nj^. 

(C) Assume G^(k 4 ) C Nj\f. By Case 3(a)(i)(B), G'^ins) There¬ 
fore, by Proposition 17.71 G^ivs) We are in the setting of 

Figure UnKa). 

(I) Assume G^(ui) C Nj^. Rule 2 does not apply to ai, hence 
G’^( 7 r) Nj\f. We are in the setting of Figure [TSl bl. The desired 
transversality follows from Propositions 15.1115.2115.31 and 15.41 
(II) Assume G^(ui) Nj\f. Assume A 4 is not of type A. Thus there 
exists u E V(fi; 4 ) n n T. By Proposition 16.81 there exists u E 
Vei^i) nAA°. By assumption u E V(k 4 ) \ {^ 1 ,^ 3 } 

We leave it to the reader to check that if A 4 in not of type A, 
then by Proposition 16.91 either vi or vs is in a contradiction. 
Thus, if G^(k 4 ) C Nj^f^ then K 4 is of type A and we are in the 
setting of Figure I15r cl and the desired transversality follows from 
Propositions 15.11 and 15.21 

(ii) Assume vs E J\f^ and vs is transparent east. Figure He) provides an 
upper bound on the possible cell types. Observe that K 3 cannot be 
of type NW or N and K 4 cannot be of type NW. Furthermore, by 
Proposition 16.91 A 3 cannot be of type NE and since the edge V 3 vi 
cannot exist A 4 cannot be of type N or NE. See Figure fTHl a). By Rule 
1 G^(a 3 ) U G^(a 4 ) U G^ivs) C Nj^. See Figure fTHlb). The existence of 
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Figure 16. (a) Possible cell types in Case 3(a) (ii). (b) Possible cell types 
in Case 3(a) (iii). (c) Possible cell types in Case 3(b). 


G^{tt) is determined by Rule 2, hence a necessary condition for G^{tt) C 
Nj^ is the existence of G^(t>i) C Nj^ or G^{v 2 ) C Nj^. 

(A) Assume G^{vi) Nj^ and G^{v 2 ) 'f- Nj^. The desired transversal- 
ity follows from Propositions 15.11 and 15.21 

(B) Assume G^{v 2 ) C Nj\f. Observe that v G Va(K 2 ). Thus Rule 2 
applied to K 2 implies that G''^(vr) C N_a. Rule 5 applies and hence 
G^(ui) C Nj^f. There are no interior boundary edges to check. 

(C) Assume G^{v 2 ) o,nd G^(r;i) C Nf/. From Figure fTHT al we 

deduce that Rule 2 does not apply to ki or K 4 . Thus G°(' 7 r) 

This puts us into the setting shown in Figure fTBT cl and the desired 
transversality follows from Propositions 15.1115.21 and 15.31 
(iii) Assume and is transparent west. Figure Etc) provides an 

upper bound on the possible cell types. Observe that K 3 cannot be of 
type NE, E, or SE and K 4 must be of type NW. Thus, we are in the 
setting of Figure Eta). By Proposition 13.41 the cell types of K 2 and K 3 
are not compatible, a contradiction. Thus this case cannot occur. 

(b) Assume V 2 G and V 2 is transparent north. Figure fTTH bl provides an upper 
bound on the possible cell types. The assumption that V 2 G AA^ forces K 2 
to be of type NE or E, and to be of type NE., N or NE. By Rule 1, 
G'^{kA) U G^{v 2 ) C N^. By Rule 2 applied to K 2 , G'^(' 7 r) C Nj^j-. Thus we are 
in the setting of Figure Etc). 

(i) Assume ^ Af^. This implies that K 3 is of type A, W, S, or SW and 
hence ^ 2,^3 G Ve(K 3 ). Since G‘^{k^) U G^{v2) U G^{i,j) C N_^f, Rule 3 
applies and hence C"'(k 3 ,^ 3 ,tt) C Nj^. Similarly, G"'(ki, ui, tt) C Nj^. 
By Proposition 13.41 ka is of type NW. This implies that vi,V 3 G Va(/« 4 ) 
Thus Rule 4 does not apply to K 4 . FigureEtc) indicates Nj^ n E{i, j)- 
That the vector field is transverse in at all the interior boundary edges 
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Figure 17. (a) Possible cell types in Case 3(a)(iii). (b) Possible cell types 
in Case 3(b). (c) Necessary grid elements in Nj^ in Case 3(b). 


(beginning with the top right and moving down to the bottom left) 
follows from Propositions 15.1115.5115.4115.4115.51 and 15.11 

(ii) Assume V3 G and V3 is transparent east. Observe that rotating clock¬ 
wise by 90° results in Case 2(b)(i) for which the desired transversality 
has been demonstrated. 

(iii) Assume V3 G and V3 is transparent west. Observe that rotating 
clockwise by 90° results in Case 2(c)(i), which has already been shown 
not to occur. 

(c) Assume V 2 G and V 2 is transparent south. We begin by rotating clockwise 
by 90° and the performing a reflection in the east-west direction. Figure [ 6 Ka) 
provides an upper bound on the possible cell types. Since V 2 is transparent 
south, k ,2 is of type SE and K 3 cannot be of type NW, N, or NE. The 
assumption that U 3 0 Af^ implies that K 3 cannot be of type E or SE This 
results in cell types as shown in Figure fTHT ai. Rule 1 implies that G^(a 3 ) U 
G^{v 2 ) C Nj^. Rule 2 implies that C Nj^. Thus the necessary set of 

tiles in Nj^ is shown in Figure [TST bi. 

(i) Assume vi G Af^ and vi is transparent north. This is the same as Case 
l(b)(i), hence the desired transversality has been demonstrated. 

(ii) Assume vi G Af^ and vi is transparent south. This is the same as Case 
2 (b) (i), hence the desired transversality has been demonstrated. 

(iii) Assume vi 0 Af^. Figure fTHf al provides an upper bound on the possible 
cell types. Since vi 0 Af^, ki cannot be of type S 01 SE. Since K 2 is of 
type SE, Proposition Id.df iil and (ix) implies that rjj G and G 
respectively. With these restrictions Proposition I3.4f i) implies that K 4 
is of type NW, N, or NE, while Proposition I3.4f xl implies that K 4 is 
of type NW, W, or SW. Thus K 4 is of type NW as is indicated in 
Figure [THlf cb 

□ 
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Figure 18. (a) The possible cell types for Case 3(c). (b) Necessary tiles 
for in Case 3(c). (c) The possible cell types for Case 3(c)(hi). 
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Figure 19. Tiles and chips for Case 3(c)(iii) 

8. Proof of Theorem 

The focus of this section is on the proof of Theorem 18.21 the main theorem of this 
paper. We begin with a preliminary result concerning the existence of global attractors 
and conclude with minor extensions of Theorem 18.21 

Proposition 8.1. Consider 

(44) X =—Tx-\-f{x), xG(0,oo)" 

where f: ( 0 , oo)"' —>■ ( 0 , oo)'^ is Lipschitz eontinuous and satisfies the conditions that there 
exist positive constants af such that 

0 < a~ < ffix) < af, i = 1 ,... ,n 

for all X G ( 0 , oo)"'. Furthermore, assume that T is a diagonal matrix with diagonal elements 
7i > 0 . Then, there exists a global compact attractor X C ( 0 , oo)”' for the flow generated 
by (j44)l . 
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We leave the proof of ProDosition IS.ll to the reader noting that for positive a sufficiently 
small and sufficiently large, the vector field (|44ll is transverse in on the boundary of 
[a-,a+]”. 

Theorem 8.2. Let = S(r, A, H^, H^) he a switching system as defined by Defi¬ 

nition \3.1\. Let J- be the associated state transition diagram given by Definition \6.3X Let 
MG(J^) = {M.{p) I p G (P, <p)} be the associated Morse graph. Choose 0 < S < 6*, where 
5* satisfies (IMl), and let 

(45) X =—Fx xG(0,oo)^ 

be an associated 6-constrained continuous switching system as defined in Section\^ Let ip 
he the flow associated with (H5D and let X C (0, oo)^ be the associated global attractor. Let 
(P, <p) be the poset given by Defi,nition \2.8[ Then, there exists a Morse decomposition for 
X under p with Morse sets M{p), p G P, and an admissible order <p. 

As indicated in the introduction Theorem 18.21 is obtained as an application of Theo¬ 
rem [2]9j Thus, the main task of this section is to establish that desired finite distributive 
lattice (N,A,V,0, 1) that satisfies hypothesis (l)-(3) of Theorem 12.91 With this in mind 
we introduce the following lemmas and proposition. 

Lemma 8.3. Let A G Att(JG). Let E{i,j) be an elementary domain with the local vector 
field in Case 2(b)(i) of the proof of Proposition [7^ If v G A, then {v,vi,V 2 } C A. 

Proof. By Figure fT^ ah ki is of type S or SE. Thus v ^ vi and hence vi G A. Further¬ 
more, K 2 is of type SE and hence by Proposition 16.91 V 2 G A. □ 

We now extend Lemma 18.81 to a sequence of elementary domains. For this need to be 
able to consider symmetric versions of Case 2(b) (i) as indicated in Figure [20l Recall that 

given an elementary domain E{i,j), |u-j, j C £{i,j). 

Lemma 8.4. Fix L > 1. Consider elementary domains E{i — i,j), i = 0,..., L. Assume 
that in region regions E{i — i,j), i = 1,..., L — 1, the possible cell types are of the form 
indicated in Figure [2W c). Assume that in region E{i — L,j) the possible cell types are 
of the form indicated in Figure \20\fa). Assume that in region E{i,j) the possible cell 
types are of the form associated with Case 3(a)(ii)(B) of the proof of Proposition Let 

D = I ^ = 0,..., l| . If G Ao, then C c Aq. 

Proof. By Lemma ISTSl G .Aq. From Figure [20llcL ^ Ve(«^ 4 ) and 

j G for K 4 associated with E{i — L -\- l,j). Similarly, v-_^jzr[ G and 

G Vai^s) for K 3 associated with E{i — L + l,j). Thus by the definition of E, 
Vi_Lj ^ and -G Therefore, j, G A- The 

result now follows by induction. □ 

Proposition 8.5. Let S = S(F, A, H^, H^) and E be as in Theorem The set 

N = {Nj^ I A G Att(A} , 
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Figure 20. (a) Possible cell types for Case 2(b)(i)* obtained from Case 
2(b) (i) by counterclockwise rotation of 90°. (b) Possible cell types for Case 
2(b)(i)** obtained from Case 2(b)(i)* by a vertical reflection, (c) Possible 
cell types in elementary domain E(i,j) in Lemma 18.41 

where is given by Definition \7.^ is a finite distributed lattice with lattiee operations 

(46) Nao a Na, ■■= NaoaAi C Nao C Na, 

(47) Nao V Nai ■■= NaoVAi = ^Ao U 
with minimal element 0 = 0 . 

Proof. Observe that N is well defined and (|46p and (1471) are well defined lattice operations 
since Att(J^) is a lattice with V and A defined by d?]) and ([ 8 ]), respectively. 

By Lemma 17.31 0 = 0. Thus all that remains to be shown is that the inclusion and 
equality of (HOl) and (|17|) hold. 

To prove that the inclusion of (1461) we consider several representative cases leaving the 
rest to the reader. Let A' := A .Ai C A.o n A.i. 

Assume that a wide chip is contained in Na'- Without loss of generality we can 
assume that this implies that there is an elementary region E{i,j) such that Na' riE{i,j) 
gives rise to Case 2(a)(i)(B) in the proof of Proposition 17.41 Observe that U 2 G W and 
hence V 2 0 A.o U A.i. If ^ AqL) Ai, then belongs to both Naq and Nai ■ If U3 € Aq 
or U 3 € Ai, then G^(u 3 ) belongs to Naq or Nai- Thus in all these cases G^(u 3 ) belongs to 
Nao^Na,. 

Assume that a narrow chip C” is contained in Na'- Again, without loss of generality 
we can assume that this implies that there is an elementary region E(i,j) such that Na' C 
E{i,j) gives rise to Case 2(a)(i)(A) or Case 3(c)(iii). Recall that in Case 2(a)(i)(A), K 3 is 
of type NE. Thus, if V 2 G Ai or V 3 G Ai, i = 0 or 1, then as indicated in Figure [S] by Rule 
1 G^(k 3 ) U G^(u 3 ) C AI 4 .. Therefore, C"' is contained in Naq 0 Nai ■ A similar argument 
applies to Case 3(c)(hi). 

Assume that a 0-tile G^{'k) is contained in Na'- If tt = (^j, rjj), then Rule 2 is applicable 
to E{i,j). Reviewing the proof of Proposition[73]we see that this happens in Cases l(a)(ii), 
1(b), l(c)(i), 2(a), 2(b), 3(b), and 3(c). In the subcases of Cases 1 and 2, v,vi G Tn A' 
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and thus u,i;i E Tn Ai, z = 0, 1. Therefore, after applying Rule 1 in these cases to iV^., 
z = 0,1, Rule 2 is applicable, and hence G°(7r) is contained in z = 0,1. The subcases 
of Case 3 follows from a similar argument based on u, z ;2 E T H 

Assume that a 1-tile G^{va) is contained in This implies that Rule 1 or Rule 

5 applies. Observe that G^{va) is introduced by Rule 1 if Ua E T n A' n T(z, j). But 
Va & TCl A' n T(z, j) implies that Ua E T 0 (A.o H Ai) ri£{i,j) and hence by Rule 1 G^{va) 
belongs to z = 0,1. If G^{vo,) is introduced by Rule 5, then without loss of generality 
we can assume that we are in the setting of Case 2(b)(i) or Case 3(a)(ii)(B). We leave it 
to the reader to check these cases using similar arguments as above. 

Finally, assume that a 2-tile G^(Ka) is contained in N_^i. This implies that Rule 0 
and/or Rule 1 applies. In both cases the fact that A1 C A.o H Ai implies that Rule 0 
and/or Rule 1 also applies to iV_ 4 ., z = 0,1. 

We now turn to the proof of equality in (jlTll . In general, since A' := AqV Ai = A-qU Ai, 
the same arguments that were used to verify the inclusion of (|46p lead to equality in the 
setting of The exception is the introduction of 1-tile G^{va) due to Rule 5. This 

occurs in two cases Case 2(b) (i) or Case 3(a)(ii)(B). 

Consider Case 2(b)(i). Since A! = AqUAi and v E A', we can without loss of generality 
assume that v E Aq- Applying Lemma IHAl we see that v,vi,V 2 E Aq and hence Rule 5 
applies to Aq. Hence A! H E{i,j) = (A.o U Ai) H E{i,j)- 

Finally, consider Case 3(a)(ii)(B). The assumption is that G^{v 2 ) C N_a', but V 2 0 Al. 
Thus, G^{v 2 ) was introduced by Rule 5 applied to E{i — 1, j). If E{i — 1, j) is of the form 
Case 3(a)(ii)(B), then we repeat the argument to conclude that G^{v 2 ) for E{i — 1, j) was 
introduced by Rule 5 applied to E{i — 2,j)- Since there are only a finite number of cases we 
can assume that there exists L > 1 such that E{i — L,j) is of the form Case 2(b)(i)* or Case 
2(b)(i)** (see Figure [20lai and (b)) and E{i — £,j) is of the form indicated in Figure ISOlci 
for .^ = 1,...,L — 1. The desired equality in (j47h now follows from Lemma 18.41 □ 

Proof of Theorem \8.2[ ProDosition l7.4| guarantees that for all A E Att(T'), E Invset^(</?). 
Thus hypothesis (1) of Theorem 12.91 is satisfied. Lemma 17^ indicates that hypothesis (2) 
of Theorem 12.91 is satisfied. Finally, Proposition 18.51 guarantees the validity of hypothesis 
(3). □ 

Remark 8.6. Because Theorem 18.21 is obtained as an application of Theorem 12.91 we also 
have information concerning the location of the Morse sets M{p) of the 5-constrained 
continuous switching system. In particular, each p E P is associated to a unique a join 
irreducible element of the lattice N and hence a unique A E Att(J-') which is join 
irreducible in Att(J-'). Theorem 12.91 guarantees that 

M{p) cNa\N^C {0,oof. 


We conclude this section by noting that the assumption that f^^^ be constant on the 
tiles is not nec essary. Applying the classical continuation theorem for Morse decompositions 
Conley 1978l | to Theorem 18.21 gives rise to the following result. 
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Theorem 8.7. Let 

( 48 ) X = —Fx + xG( 0 ,oo)^ 

be a 5-constrained continuous switching system with an associated Morse graph MG for the 
global attractor of ( 1481 ) . If 

sup ||/(x) -/(‘^)(x)|| < e 

a;g(0,oo)2 

for sufficiently small e > 0, then MG zs a Morse graph for the global attractor of 

x = —Fx + f{x), xG(0,oo)^. 


9. Conclusion 


We conclude with a few comments concerning the results of this paper and potential 
future directions. 

We begin by reviewing the conclusion of Theorem 18.21 and in particular the Morse set 
M (p). Extending the discussion of Remark 18.61 note that 

1= V Mg N 

.AGAtt{J') 

and hence 

M{p) C X\iVi. 

Whether or not M{p) 0 is depends on the system being considered. However, given 
the construction of T this cannot be determined given the calculations presented in this 
paper. This is essentially due to the fact that we do not identify any dynamics that is 
associated with white vertices. We are currently working on algorithms for generating an 
alternative state transition diagrams that are capable of capturing the dynamics imposed on 
(5-constrained continuous switching system by white vertices. However, the computational 
cost of working with these alternative state transition diagrams is greater. Whether the 
additional dynamical information that can be gained is worth the additional cost in the 
context of the analysis of biologically motivated networks remains to be seen. 

With regard to (5-constrained continuous switching systems there are three points worth 
emphasizing. First, within the b collars we impose minimal assumptions on the form of 
Thus we do not need to assume that f^^^ is based on a particular nonlinearity, e.g. 
a Hill function, and thus the dynamics we are recovering is valid for an extremely wide 
set of potential models. Second, the larger b is the less steep f^^^ need be, and hence, 
the less switch-like the system needs to be. Third, we give explicit bounds on b* in terms 
of the parameters of A and T. A possible consequence of this is that the computational 
tools developed for switching systems can be used to guide the study of the local and 
global dynamics of arbitrary systems of the form ([T]) for fixed families of nonlinearities. 
To be more precise, assume that the /„ are given in terms of Hill functions and we are 
interested in particular dynamical structures. Letting the Hill coefficients /c —)• oo produces 
a switching sy stem, the global dynamics of which can be analyzed over all of parameter 
space using T Cummins et al.l |2ni5b |. We then can identify parameter values at which 
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the desired nonlinear dynamics is exhibit and determine the maximal size of perturbation 
6 *. Given 5* one can choose k sufficiently large so that the Hill function approximates a 
(5-constrained continuous switching system for 6 < 6*. More standard numerical methods 
can then be used to identify the desired Morse set for ([1]) with this large Hill coefficient k. 
Finally, numerical continuation techniques can be employed to determine if the dynamics 
continues to lower biologically motivated values of k. 

The focus of this paper is on translating information obtained from piecewise constant 
models in the form of switching systems that are motivated by regulatory networks, to 
information about the dynamics generated by Lipschitz continuous differential equations. 
However, as is mentioned in the introduction part of the motivation for this paper is 
our interest in the mathematically rigorous analysis of global dynamics for multiparameter 
systems over large regions of parameter space and the question of whether for these purposes 
it is computationally efficient to use the techniques described here. With this in mind the 
results presented in Section [7] are much more general than those required to study regulatory 
networks with two nodes. Observe that Theorem 18.71 suggests that given a systems of 
ordinary differential equations of the form ([T]) one could try to compute the associated 
dynamics by approximating / via a linear term T and a piecewise constant function A, 
and then, computing the associated Morse graph for the associated switching system. We 
plan to exploring the effectiveness of such a procedure. However, there are at least two 
related issues that need to be addressed. First, we need explicit results for bounds on e, 
the acceptable size of perturbation in Theorem 18.71 Second, we need to understand how 
to determine the threshold values used to defined the domains of the piecewise constant 
functions. 
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